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SECTION  I 
INTRODUCTION 


The  understanding  and  modelling  of  plasma  convection 
in  the  high  altitude  environment  is  a  central  issue  in  the 
proper  simulation  of  a  realistic  nuclear  scenario.  This 
paper  is  principally  concerned  with  the  study  of  a  key  aspect 
of  the  problem  which  we  have  chosen  to  label  "striation 
convection".  One  distinguishes  between  usual  or  macroscopic 
plasma  convection  and  striation  convection  because  of  present 
day  limitations  in  computational  techniques  and  resources. 

In  principle,  there  is  only  a  single  physics  of  plasma 
motion  governed  by  the  same  hydrodynamics,  whether  the 
scale  is  1000  km  or  1000  cm.  In  practice,  however,  we  are 
forced  to  arbitrarily  divide  the  phenomena  into  manageable 
scale  size  regimes.  We  treat  the  large  dimensions  (limited 
by  the  computational  grid)  directly  with  a  suitable  plasma 
fluid  code.  For  late-time  studies,  an  electrostatic  routine 
is  the  usual  choice.  Shorter  scales  are  defined  as 
"striations"  and  are  modelled  in  various  types  of  auxiliary 
routines.  In  previous  research  and  systems  analysis,  there 
has  been  no  provision  for  an  independent  motion  or  convection 
of  this  small  scale  plasma.  While  some  previous  work  has 
ignored  the  motion  of  this  plasma  altogether,  even  the  most 
advanced  computations  have  only  considered  that  it  will  move 
in  an  identical  manner  to  the  smallest  scale  plasma  that  the 
grid  does,  in  fact,  compute.  That  is,  it  is  assumed  that  the 
striations  are  "imbedded"  in  the  local  gross  scale  plasma  and 
will  be  carried  along  by  that  plasma  in  a  Lagrangian  manner. 

Even  if  one  were  unaware  of  the  detailed  physics  of 
the  gradient  drift  mechanism,  the  present  computations  would 
appear  to  contain  an  inconsistency.  We  would  observe  from  the 
numerical  data  that  plasma  features  of  different  densities 


would  be  convecting  at  different  rates.  It  would  thus  seem 
intuitively  incorrect  that  all  features  below  the  grid 
dimension  should  be  treated  as  a  common  Lagrangian  mass 
with  uniform  motion. 

If  we  turn  to  a  direct  consideration  of  the  gradient 
drift  physics,  it  is  immediately  apparent  that  to  neglect  the 
relative  convection  of  small  scale  plasma  is  to  actually 
neglect  the  gradient  drift  mechanism  itself.  It  is,  after 
all,  the  relative  slip  of  different  features  that  produces 
structure  on  the  fine  scale  in  the  first  place. 

At  late  times,  the  process  of  striation  convection 
has  profound  effects  upon  the  overall  plasma  evolution. 

Both  the  breakup  of  large  initial  features  and  the  dispersal 
of  material  to  remote  locations  is  governed  by  the  mechanism. 
It  could  easily  be  argued  that  the  process  is  central  to 
all  late-time  modelling  of  plasma  evolution.  Clearly, 
without  a  proper  consideration  of  this  convection,  all  con¬ 
ventional  computations  of  high-altitude  plasma  motion  lead 
to  unphysical  results  at  relatively  early  times. 

In  lieu  of  the  ability  to  carry  out  "exact"  plasma 
simulations  over  the  vast  range  of  scale  sizes  actually 
encompassed  in  a  nuclear  scenario,  our  proposal  is  to  model 
the  plasma  in  an  integral  manner  by  conserving  certain  key 
extensive  properties.  In  particular,  the  modelling  will  rely 
on  conservation  of  plasma  density  and  plasma  "power".  Power 
is  defined  as  a  function  of  the  square  of  .the  density,  either 
absolute  or  as  a  variance  about  a  local  mean.  The  next 
section  of  the  paper  will  explicitly  review  the  construction 
of  the  class  of  power  functions  to  be  employed.  The  use  of 
plasma  power  to  model  convection  is  closely  tied  to  previous 
work  on  power  modelling  of  in-situ  microstructure.  We  will 
reference^  and  make  use  of  previous  results  in  constructing 
the  present  model. 
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Fundamentally,  the  physics  of  striation  convection 
produces  a  differential  slip  of  high  and  low  density  plasma 
within  a  macroscopic  cell.  The  result  is  a  convective  inter¬ 
leaving  of  high  and  low  density  which  appears  as  "mixing" 
on  the  grand  scale.  Mixing  is,  of  course,  a  diffusive 
process  and,  thus,  the  apparent  end  result  is  displayed  as 
plasma  diffusion  among  neighboring  macroscopic  cells.  The 
process  is  not  to  be  confused  with  true  crossfield  diffusion, 
but  rather  is  the  integrated  result  of  plasma  striation  inter¬ 
leaving. 

The  next  section  of  the  paper  will  develop  the  concept 
of  plasma  power  and  lay  the  groundwork  for  the  detailed 
convection  model.  A  subsequent  section  will  outline  the 
actual  mixing  model  and  indicate  the  connection  to  previously 
developed  microstructure  modelling.  An  appendix  will  display 
a  simple  microstructure  model  which  incorporates  recent 
philosophical  changes  in  the  older  model  and  provides  a 
convenient  vehicle  for  testing  the  mixing  concept.  The  final 
sections  give  example  computations  and  conclusions. 
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SECTION  II 

PLASMA  POWER  FORMULATION 

It  is  well  known  that  field  line  integrated  plasma 
density  in  the  F-layer  of  the  ionosphere  is  governed  by 

2 

incompressible  fluid  dynamics  for  crossfield  convection. 

The  relative  uniformity  of  plasma  parallel  to  the  field 
allows  a  reasonable  two-dimensional  approximation  for 
modelling  transverse  microstructure  phenomenology.  Thence, 
with  the  exception  of  diffusion  on  the  very  smallest  scale 
size,  the  plasma  density  itself  can  be  modelled  as  incom¬ 
pressible  flow.  It  follows  that,  in  addition  to  the  density, 
any  function  of  density  becomes  a  preserved  Lagrangian 
quantity.  Of  particular  interest  is  the  square  of  plasma 
density  which  is  equivalent  to  "plasma  power".  Previous 
research1  has  laid  the  groundwork  for  the  construction  of 
the  conservation  concepts  to  be  employed  here. 

For  convenience,  let  us  consider  a  two-dimensional 
Cartesian  space  containing  a  plasma  density  distribution. 

The  space  is  sufficiently  large  that  the  plasma  density  is 
uniform  at  the  boundary  for  all  time  periods  of  interest. 

The  schematic  in  Figure  1  displays  an  actual  continuous 
density,  n(x,y)  ,  by  the  dashed  line  and  a  finite  cell 
representation,  nQ(x,y)  .  Denoting  the  cell  dimension  as 
2L  and  the  total  spatial  extent  as  ML,  we  can  construct  a 
variety  of  conservation  relations.  For  the  density  itself, 
we  can  express  the  continuous  function  as 
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Figure  1.  Background  power. 
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In  thermodynamic  terminology,  we  shall  find  that  density 
conservation  in  gradient  drift  modelling  is  analogous  to 
"energy"  in  conventional  fluid  turbulence  modelling. 

For  the  density  square  relations,  we  can  write 


JL 
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-total  diffusion  1  0 


local  diffusion  (±) 
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To  continue  the  thermodynamic  analogy,  we  will  find  that 
density  square  conservation  resembles  "inverse  entropy". 

To  introduce  the  concept  of  plasma  power,  let  us 
represent  the  density  as  a  finite  Fourier  series 
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where  a  is  the  real  polar  modulus  and  0  is  the  plasma 
phase.  The  density  square  representation  becomes 


M2L 


r~2  jjn2dx  dy  =  ^  a2(kx,ky)  =  EPSD 
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x '  y 


Thence,  the  conservation  relation  which  governs  all 
microstructure  evolution  becomes 


D  I  local  j 

Dt  PSD^2L  "  [diffusion] 


2L 


(7) 


where  EPSD  is,  of  course,  the  power  spectral  summation. 
When  referenced  to  a  local  mean  or  detrending  curve,  it  is 
commonly  called  the  density  variance. 
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An  exact,  continuous  treatment  of  the  evolution  of 
plasma  power  would  provide  no  simplification  over  an  exact 
density  solution  for  convection.  It  is  the  ability  to 
divide  and  approximately  model  various  components  of  the 
power  which  makes  the  general  concept  useful. ^ 

The  basic  component  of  power  is  provided  by  the 
macroscopic  background.  Again,  with  reference  to  Figure  1, 
we  can  define  both  a  local  background  power,  ,  and  a 

total  background  power,  P 


n^(x,y)  4L2 


(8) 


M 


■I 


(9) 


Note  that,  in  the  limit  of  infinitely  small  cell 


goes  to  zero. 


size,  2L  ,  the  absolute  value  of  I 
2 

However,  M  becomes  proportionately  infinitely  large  and 
PT  would  give  the  true  total  power  in  the  overall  space. 

For  finite  L  ,  the  expression  PT  gives  only  the  "back¬ 
ground  power"  with  several  components  remaining  to  give 
the  true  total. 

In  Figure  2,  we  focus  on  an  individual  cell  in  order 
to  define  components  of  the  "local  power".  On  the  left, 
we  describe  the  "local  macrb  power",  P^  ,  using  a  linear 
detrending  line  of  slope,  d  .  Note  that  the  line  is  drawn 
to  provide  no  net  integrated  density.  This  construction  is 
important  to  avoid  cross  terms  in  defining  the  power  and  ease 
in  later  adding  components  of  power  together. 
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is  the  local  slope  power  defined  as 


where  F 

o 


F 


o 


(ID 


Aside  from  its  role  in  detrending,  the  quantity  Fq 
is  important  because  we  identify  it  as  the  local  driver 
of  plasma  microstructure  within  the  cell.  To  make  a  connection 
with  previous  work,1  we  can  rewrite  (11)  in  terms  of  a  grid 
cell  reference  wavenumber  k  as 


(12) 


Note  that  kQ/Fo  is  independent  of  cell  size,  as  one  would 
expect  of  a  true  physical  quantity.  The  value  of  Fq  , 
however,  will  disappear  and  play  no  role  in  power  conservation, 
if  L  is  made  infinitely  small.  Thus,  Fq  is  a  model- 
dependent  artifice  which  provides  one  component  of  the  density 
variance  within  a  grid  cell. 

For  illustration  purposes,  we  have  oriented  the  slope 
along  a  particular  face  of  a  square  cross-section  cell. 
Obviously,  in  general,  a  given  cell  will  demand  a  detrending 
slope  surface  at  some  angle  to  the  natural  cross  section.  In 
application,  the  modelling  allows  an  arbitrary  orientation  of 
the  slope  surface.  The  peculiar  intersections  between  cells 
are  not  considered,  of  course,  and  provide  one  of  several 
sources  of  inconsistency  common  to  finite  difference  formu¬ 
lations  . 

On  the  right-hand  side  of  Figure  2,  we  focus  on  the 
remaining  plasma  power  defined  as  "local  micro  power".  This 
quantity  is  treated  as  a  normalized  power  spectral  summation 
denoted  by  Fc  . 
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(13) 


pc  =  I  f2 

The  general  formulation  and  spectral  characteristic  of 
this  quantity  have  been  discussed  in  previous  work.^  The 
specific  simplifications  that  will  be  employed  in  the 
present  study  are  shown  in  the  Appendix.  The  major  assumption 
concerning  this  component  of  power  is  the  use  of  a  single 
one-dimensional  spectrum  with  modulus  amplitudes,  f  ,  to 
represent  all  orientations  of  the  general  transform  with 
modulus  amplitudes,  a  . 

In  summary,  we  note  that  the  total  power  in  a  cell, 

P  ,  is  given  by 

I 

P£  =  P4(l  +  Fq  +  Fc)  (14) 

and  the  total  power  in  the  entire  space,  PT  ,  is  given  by 


Again,  it  is  important  to  remember  that  the  components 
Fq  and  Fc  arise  only  because  of  finite  gridding. 

In  the  next  section,  we  will  employ  this  plasma  power 
formulation  to  construct  a  microstructure  convection  model. 
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SECTION  III 
THE  MIXING  MODEL 


The  gradient  drift  convection  of  microscopic  features 
produces  an  effective  interleaving  of  plasma  between  neigh¬ 
boring  macroscopic  cells,  as  shown  in  Figure  3.  When  viewed 
from  the  standpoint  of  macroscopic  plasma  power,  the  net 
effect  of  the  interleaving  is  to  reduce  the  combined  power 
of  the  cells.  The  simple  schematic  of  Figure  4  illustrates 
the  point. 

In  part  A  of  Figure  4,  we  display  a  simple  pair  of 
adjacent  cells  before  convection.  Each  has  the  same  area,  A  , 
but  different  and  uniform  background  densities,  n^  and  n2  . 
In  part  B,  we  interchange  equal  pieces  of  area  A'  and 
corresponding  densities,  commensurate  with  imcompressible 
microscopic  interleaving.  The  proper  cell  boundary  is  now 
structured  as  shown.  Finally,  in  part  C,  we  redraw  the 
original  grid  boundary  and  construct  new  average  or  "mixed" 
densities,  n^  and  n^  • 

Prior  to  interleaving  the  combined  cell  power,  Pc 
may  be  stated  as 


After  interleaving,  the  power  becomes 

Pi  =(ni2  +  n22)4t2 

=  (n-^  -  6n1  +  6n2)2  4l2 


where  6 


+  (n2  -  6n2  +  6n1) 2  4L2 
A '/A  and  A  =  4L2  . 


(16) 


(17) 
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The  change  in  macroscopic  power  APC 

APc  =  P£  -  Pc  =  -2(6-62) (n1-n2)2  4L2  (18) 

We  should  note  that  geometry  prevents  6  from  exceeding 
1/2  (the  maximum  power  loss  in  the  pair) . 

Initially,  the  defined  microscopic  power,  PM  ,  is 
zero.  After  mixing,  referenced  to  the  new  primed  macroscopic 
cell  densities,  it  would  be  computed  as 

(n_-n* ) 2  A'  +  (n.-n*)2  (A-A')  , 

P'  -  -2-1 - _L_1 -  4L2 

(n. -ni) 2  A '  +  (n--nl)2  (A-A*)  ~ 

+  — - - - - ^ - - - - -  (19) 

Substitution  and  reordering  reduces  to 

P^  =  2(6-62) (ni-n2)2  4L2  (20) 

The  change  in  microscopic  power  is,  thus,  identically 

APM  -  PA  -  PM  -  _APC  (21) 

Overall  total  power  has  been  conserved,  but  macroscopic  power 
has  been  lost  in  the  exact  amount  equal  to  microscopic  power 
gain.  This  equality  and  conservation  forms  the  basis  for  the 
mixing  model. 

We  have  constructed  a  simple,  but  physically  reasonable, 
mixing  model  by  considering  the  interaction  between  a  cell 
and  its  closest  four  neighbors.  The  five-point  "cross"  that 
results  is  shown  in  Figure  5.  Eventually,  it  may  be  necessary 
to  employ  a  more  elaborate  13-point  model  to  provide  better 
symmetry  and  properly  treat  the  slope  power  component.  For 


Figure  5.  Simple  mixing  model. 


the  present,  however,  we  believe  that  the  simple  model  is 
adequate  to  account  for  the  basic  physics  and  explore  the 
resultant  phenomenology. 

The  total  background  power  in  the  cells  of  Figure  5 
is  given  by 


P  =  4L2 (n?  .  +  n2  .  +  n2  .  .  +  n?  .  +  n2  .  ) 

C  1,3  1+1,3  i“l»3  1,3+1  1,3-1 


(22) 


Prior  to  a  mixing  step,  we  can  define  the  density  differences 
between  cells  as 


An  =  (n 


i/3 
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(23) 


After  a  mixing  step,  the  new  cell  densities  may  be  expressed 
as 
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(24) 


where  the  transfer  fractions,  T  and  T  are  analogous  to 

x  y 

the  simple  fraction  6  in  the  earlier  example.  The  actual 
mixing,  of  course,  occurs  between  cells  aligned  with  the  slip 
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velocity  vector  between  neutral  and  ion  motion.  The  quantities 
and  are  merely  components  defined  for  convenience  in 

the  specific  coordinate  system.  We  can  readily  define  a  single 
scalar  mixing  parameter,  a  ,  by  relating  the  transfer  to 
the  slip  components. 


The  U's  are  the  macroscopic  plasma  velocity  as 

2 

determined  by  the  usual  macroscopic  potential  solver  and 
the  V's  are  the  neutral  wind  velocity. 

Substituting  (23)  and  (25)  into  (24)  and  thence  into 
(22)  gives  the  new  cell  power  after  a  mixing  step.  Subtracting 
this  power  from  the  initial  power  provides  an  expression  for 
the  macroscopic  power  loss  during  a  time  step. 

-APC  =  ^2S2a2  ( A2  +  A2  +  A]_A2)  +  2C2ct2(A2  +  A2  +  A^) 

+  2SCa2(A1  +  A2) (A3  +  A  )  -  2Sct(A2  +  A2) 

-  2Ca  (A2  +  A2)j4L2  (26) 

Now,  for  each  time  step  for  each  cell,  the  quantity 
-AP_(x,y)  is  a  known  distribution  from  the  solution  of  the 

,  i 

microstructure  routine.  For  the  example  calculations  to  be 
presented  in  the  next  section,  it  is  easily  expressed  from  the 
simple  model  outlined  in  the  Appendix. 
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-Apc  =  APm 
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1  o  o  o 


At 


(27) 


where  At  is  the  time  step. 

At  each  time  step,  the  relation  (26)  -  (27)  is  solved 

to  obtain  a(x,y)  and  thence  transfer  functions  T  and 

x 

T  constructed.  Both  the  density  distribution,  n  ,  and 
the  microscopic  power  distribution,  F  ,  are  thence  mixed 
to  complete  the  operation. 

A  complete  time  step  for  the  plasma  consists  of  the 
following  four  operations. 


1.  Potential  solution  to  obtain  U  ,  U 

x  y 

2.  Microscopic  power  computation 

3.  Mixing  of  density  and  microscopic  power 

4.  Macroscopic  convection 

In  the  next  section,  we  will  display  a  series  of 
sample  computations  and  discuss  the  results. 
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SECTION  IV 

EXAMPLE  CALCULATIONS  AND  RESULTS 

One  of  the  classic  problems  in  numerical  simulation 
of  gradient  drift  phenomenology  is  the  axisymmetric  plasma 
distribution  and  wind  pattern  created  by  a  fireball  in  a 
vertical  magnetic  field.  In  the  absence  of  striation  con¬ 
vection,  there  is  no  plasma  dynamic  mechanism  for  breaking 
up  and  dispersing  the  cloud.  The  example  calculations,  to  be 
presented  here,  will  focus  on  a  simple  representation  of 
that  geometry  to  display  the  behavior  of  the  convection  model 
when  faced  with  this  classic  configuration. 

Each  of  the  three  examples  will  begin  with  an  initial 
Gaussian  plasma  density  distribution  of  the  form 


where  the  density,  n  ,  is  obviously  normalized  to  a  uniform 
background  and  the  r's  are  the  various  radii.  In  all 
cases,  the  Gaussian  radius,  rQ  ,  will  be  100  units.  The 
unit  dimension  is  arbitrary,  but  may  be  thought  of  as  a  kilo¬ 
meter  for  convenience.  The  center  position,  r^  ,  will  be 
zero  for  the  first  and  third  cases,  but  offset  100  units  in 
the  second  case  . 

The  computational  grid  will  be  a  Cartesian  (x,y)  space 
of  total  extent  1000  units  on  a  side  divided  32x32.  The  grid 
reference  dimension,  L  ,  is,  thus,  15.6  units.  These  numbers 
provide  an  outer  scale  reference  wavenumber,  kR  ,  of  .22 
and  an  inner  scale  reference  wavenumber,  kD  ,  of  220.  All 
problems  commence  with  the  microstructure  power,  Fc  ,  set 
equal  to  zero. 


21 


The  wind  pattern  is  fixed  in  time  and  the  same  for 
each  case. 


V 


r 

r  +  r 
o 


(29) 


The  wind  is,  thus,  radially  symmetric  with  the  zero  value  at 
the  center  and  growing  to  one  unit  per  second  at  large  radii. 
The  macroscopic  potential  is  solved  in  the  Neumann  mode 
which  transmits  plasma  across  the  boundary  in  the  wind 
direction. 

All  time  steps  are  ten  seconds,  which  implies  three 
steps  per  2L  for  expected  peak  velocities  of  one  unit  per 
second.  The  time  interval  appears  adequate  for  the  first 
two  cases,  but  is  perhaps  a  little  coarse  for  the  third  case 
(as  the  results  will  indicate) .  The  various  mesh  reductions, 
necessary  for  each  of  the  steps  corresponding  to  potential 
solution,  mixing,  and  convection,  utilize  alternating  sweep 
techniques.  Unphysical,  numerical  asymmetries  can  creep 
into  the  pattern  during  a  coarse  time  step  in  spite  of  our 
consistent  use  of  predictor-corrector  procedures. 

In  the  first  case,  the  cloud  is  placed  on  center  with 
the  wind.  The  initial  density  contours  are  as  shown  in 
Figure  6A.  In  all  of  the  2-D  plots,  the  contour  values  are 
shown  above  the  figure  and  the  mesh  of  letters  A  to  I 
represent  intermediate  levels  between  ascending  contours. 

In  the  purely  axisymmetric  case,  there  is  very  little 
plasma  motion  until  the  microstructure  has  had  a  chance  to 
build  to  significant  values.  The  first  noticeable  effects 
occur  around  600  seconds.  Figures  6B  and  6C  display  contours 
of  n  and  Fc  .  The  toroidal  distribution  of  Fc  is  the 
expected  pattern,  which  gives  a  ring  of  peak  values  on  the 
steep  outer  face  of  the  cloud. 
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Figure  6B.  Case  I  -  plasma  density  -  600  seconds. 
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Visualization  of  the  plasma  evolution  is  easiest  using 
3-D  plots  and,  thus,  will  be  employed  for  all  of  the  later 
time  development.  Figure  6D  at  100  seconds  shows  the  density 
in  3-D  prior  to  any  significant  motion.  Figure  6E  at 
1000  seconds  reflects  the  initiation  of  noticeable  mixing. 

The  finite  cell  size  and  the  lack  of  radial  symmetry  create 
the  four  prongs  to  the  structure.  If  the  gross  plasma  had 
possessed  any  net  macroscopic  velocity  (as  in  the  next  case) , 
the  prongs  would  have  been  washed  out.  As  it  is,  they 
emphasize  the  "cross  pattern"  of  the  mixing  model.  Figures 
6F  and  6G,  at  1300  and  1500  seconds,  show  the  further 
development  of  structure.  The  decline  of  peak  density 
in  the  central  cloud  is  evident  by  this  time.  By  2000 
seconds  in  Figure  6H,  a  significant  dispersal  of  plasma 
has  occurred  and  material  is  reaching  the  grid  boundary.  The 
remaining  plots.  Figures  61,  6J,  and  6K,  at  2500,  3500, 
and  6000  seconds,  display  the  final  evolution  of  the  plasma. 
There  is  complete  destruction  of  the  original  cloud  and 
dispersal  of  plasma  beyond  the  grid  boundaries. 

Fundamentally,  the  first  case  demonstrates  the  overall 
behavior  that  one  expects  from  a  striation  convection  model. 
There  are,  however,  certain  obvious  deficiencies  with  respect 
to  the  details  of  the  computation.  The  stationary  character 
of  the  initial  plasma  and  the  lack  of  complete  symmetry  in 
the  model  combine  to  produce  the  cross  pattern  in  the  dispersal. 
Once  it  is  generated,  the  cross  pattern  appears  to  augment 
itself  by  serving  as  a  perturbation  in  the  macroscopic 
potential  solver.  In  a  sense,  the  macroscopic  fingers,  once 
produced,  may  be  outrunning  the  mixing  process.  The  secondary 
density  peaks  near  the  boundary,  as  in  Figures  6H  and  61,  would 
indicate  this  influence.  On  the  other  hand,  the  mixing  does 
not  completely  lose  its  effect,  as  evidenced  by  the  continuous 
decline  of  the  peak  central  density. 
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Figure  6F .  Case  I  -  plasma  density  -  1300  seconds. 
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Figure  61.  Case  I  -  plasma  density  -  2500  seconds. 


UBIWMU'W.  «■ 


There  appear  to  be  three  issues  raised  by  the  results 
of  the  first  case  which  require  further  thought.  First,  if 
the  mixing  model  were  absolutely  symmetric,  the  macroscopic 
fingers  would  not  have  occurred.  Thus,  future  work  needs 
to  evaluate  the  importance  of  using  the  more  elaborate 
13-point  model  mentioned  in  an  earlier  section.  On  the  other 
hand,  no  finite  difference  procedure  can  be  absolutely 
"smooth"  and,  perhaps,  the  results  are  merely  exaggerating 
a  discrepancy  in  the  modelling  which  needs  specific  attention. 
In  particular,  we  identify  a  second  topic  as  the  question  of 
whether  the  macroscopic  motion  should  be  able  to  "outrun" 
the  mixing  regardless  of  the  perturbation  produced.  Concep¬ 
tually,  it  seems  that  the  "mixing  front"  ought  to  proceed 
at  about  the  same  pace  as  a  "macroscopic  finger  front".  The 
two  are,  in  fact,  supposed  to  represent  the  same  physical 
phenomenon  and  the  specific  coupling  to  the  wind  is  not 
dependent  on  scale  size  to  first  order. 

The  third  issue  raised  by  the  results  of  the  first 
case  is  the  virtual  uniqueness  of  the  initial  condition, 
which  produces  no  gross  macroscopic  motion.  In  any  real 
scenario,  there  would  not  be  complete  axisymmetry  and  the 
gross  motion  would  undoubtedly  wash  out  perturbations  produced 
by  the  mixing  model.  The  purpose  of  case  two  is  to  explore 
this  issue  in  more  detail. 

In  the  second  case,  we  start  the  computation  with  the 
same  basic  wind  and  cloud,  but  now  the  cloud  is  offset 
diagonally  from  the  center  by  100  units.  The  initial  plasma 
distribution  is  seen  in  Figure  7A.  By  600  seconds,  a 
significant  macroscopic  velocity  (caused  by  the  offset)  has 
pushed  the  cloud  into  a  barium  release  type  configuration,  as 
seen  in  Figure  7B.  The  microstructure  power  distribution  at 
this  time,  as  displayed  in  Figure  7 C,  is  quite  different  than 
the  earlier  case.  As  expected,  however,  its  intense  region 
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is  centered  on  the  edge  pointed  with  the  neutral  wind.  The 
convection  of  microstructure,  forming  tip  wakes,  is  evident 
in  the  figure.  Figure  7D  is  the  3-D  plasma  plot  at  100 
seconds,  whereas  Figure  7E  is  the  further  evolution  at 
400  seconds. 

At  800  seconds,  in  Figure  7F,  the  cloud  has  reached 
the  boundary  and  is  starting  to  leave  the  grid.  We  should 
note  two  things  about  the  figure.  First,  as  expected,  there 
are  no  "cross-like"  perturbations.  The  macroscopic  motion 
has  clearly  washed  them  out.  Second,  at  this  last  time 
frame,  before  the  cloud  leaves  the  grid,  there  is  not  much 
evidence  of  peak  density  decay.  It  is  not  clear  that  the 
mixing  has  played  a  significant  role  in  the  cloud  evolution. 
Whether  it  should  have  played  a  role,  at  this  particular 
time,  is  not  clear  from  our  experimental  knowledge  of  such 
clouds.  We  do  wonder,  however,  as  in  case  one,  if  the 
mixing  has  been  "too  weak"  to  keep  pace  with  the  macroscopic 
motion. 

Figures  7G,  7H,  and  71  display  the  ultimate  macroscopic 
convection  of  the  case  two  cloud  out  of  the  grid.  While 
interesting  in  displaying  the  ability  of  the  numerical  tech¬ 
niques  to  handle  the  task,  they  shed  little  light  on  the 
specific  influence  of  striation  convection. 

In  the  course  of  the. present  work,  we  have  simultaneously 
begun  the  research  and  development  on  the  new  microstructure 
model  mentioned  in  the  Appendix.  One  of  the  early  observations 
was  that,  while  the  simple  model  generally  agreed  in  quali¬ 
tative  format  with  the  new  results,  it  consistently  produces 
much  lower  levels  of  microstructure.  Thus,  it  occurred  to 
us  that  we  ought  to  try  the  mixing  routine  using  artificially 
augmented  microstructure  rates  in  the  simple  model.  Using 
some  modifications  to  the  factors  6^  and  B3  >  we  can  make 
the  simple  microstructure  model  come  closer  to  the  new  research 


39 


DENSITY 


results.  (Cases  I  and  II  were  run  with  (3^  =  6^  =  !•)  For 
purposes  of  the  present  work,  the  details  are  unimportant, 
but  the  fact  of  augmentation  in  as  a  qualitative  feature 

is  the  issue.  The  third  case  is  a  rerun  of  case  one  with 
artificially  augmented  Fc  (3  =  10). 

Figure  8A  is  a  2-D  plot  of  plasma  density  in  Case  III 
at  100  seconds.  We  may  note  that  the  peak  central  density 
has  already  been  reduced  to  4.97  from  5.00.  Of  prime 
interest  is  Figure  8B,  which  shows  Fc  at  100  seconds. 
Comparison  with  Figure  6C  for  Case  I  at  600  seconds  indicates 
that  the  peak  values  for  Fc  are  more  than  seven  times 
larger,  even  though  the  time  is  six  times  shorter.  The  3-D 
density  plot  at  100  seconds  is  shown  in  Figure  8C. 

The  evolution  of  plasma  in  Case  III  is  very  rapid,  as 
displayed  in  Figures  8D,  8E,  8F,  and  8G  at  200,  300,  400, 
and  500  seconds  respectively.  The  artificially  high  rate 
of  microstructure  growth  speeds  up  the  mixing  and  rapidly 
dissipates  the  cloud.  Even  in  the  stationary  environment, 
the  process  is  relatively  smooth  compared  to  the  earlier 
Case  I.  Clearly,  the  mixing  has  completely  overwhelmed 
the  macroscopic  motion  in  this  example.  While  there  is  some 
evidence  of  the  cross  pattern  in  Case  III,  we  would  conclude 
that  the  bulk  of  the  effect  in  Case  I  was  due  to  augmentation 
by  the  potential  solver. 

The  tuning  of  future  microstructure  models  to  provide 
consistent  motion  between  the  striation  convection  front  and 
the  macroscopic  feature  front  is  probably  an  important 
consideration . 
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SECTION  V 
CONCLUSIONS 


The  use  of  power  conservation  to  construct  a  mixing 
model  representation  of  striation  convection  appears  to  be 
quite  successful.  The  qualitative  features  are  in  agreement 
with  the  expected  behavior  of  such  physics.  With  the  future 
use  of  an  advanced  microstructure  routine  and  integration 
with  experimental  data,  we  believe  that  satisfactory  quanti¬ 
tative  results  are  possible.  The  procedure  will  shortly  be 
added  to  the  SCENARIO  systems  analysis  code  and  provide  a 
great  increase  in  ability  to  run  to  late  times.  We  would 
recommend  that  the  model  should  be  added  to  research  codes 
such  as  MELT  and  MAGIC. 

As  indicated  in  the  commentary  of  the  text,  there  are 
a  variety  of  studies  and  tasks  which  would  enhance  the 
understanding  and  general  credibility  of  the  technique. 

At  least  for  research  purposes,  the  development  and  utili¬ 
zation  of  a  more  symmetric  13-point  model  ought  to  be  explored. 
The  chief  area  of  work,  however,  should  focus  on  the  close 
dependence  with  microstructure  rates  and  the  implications  in 
more  general  scenarios. 

Clearly,  the  detailed  research  and  development  required 
for  deployment  in  the  full  three-dimensional  environment 
remains  relatively  unexplored.  Striation  convection  under 
the  influence  of  gravity  and  inertial  forces  are  topics  for 
future  work. 
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SECTION  VII 

APPENDIX  -  SIMPLE  MICROSTRUCTURE  MODEL 

Initially,  when  the  construction  of  the  mixing  model 
was  undertaken,  it  was  assumed  that  it  would  be  coupled 
directly  to  the  microstructure  model  of  Reference  1.  It 
appeared  reasonable,  however,  to  make  available  a  simple 
model  for  both  quick  running  and  easily  visualized  behavior 
during  code  construction  and  refinement.  Recently,  however, 
guided  by  some  astute  observations  on  the  true  physics  of 
outer  scale  behavior  by  Wittwer ,  we  have  undertaken  a 
significant  revision  of  the  initial  code.  Interestingly 
enough,  the  simple  model  to  be  described  here  is  phenomeno¬ 
logically  closer  to  the  new  code  than  the  old.  The  end 
result  of  all  of  these  considerations  is  that  we  have  chosen 
to  use  the  simple  code  in  generating  the  results  in  the 
present  study  instead  of  the  more  complex  model  of  Reference  1. 

Assuming  that  the  reader  is  either  familiar  with  the 
earlier  work  or  can  make  reference  to  it , ^  the  construction 
of  the  simple  model  is  straightforward.  Essentially,  we 
assume  that  both  the  outer  scale  and  inner  scale  are  fixed. 
Thence,  we  consider  only  the  cascade  regime  of  power.  That 
is,  all  of  the  microstructure  power  to  be  treated  in  the 
problem  is  assumed  to  be  .  There  is  only  one  rate  relation 

and  it  refers  to  this  specific  regime  (Equation  79  of 
Reference  1) . 
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The  constitutive  relations  then  become 


*C  =  /kDkR 

kD  =  (5/3) %  kD  (A-2) 

kD  *  QkR 

where  kR  is  now  the  fixed  outer  scale  wavenumber,  kQ  , 
corresponding  to  the  cell  grid  dimension  parameter 


The  inner  scale  wavenumber,  kQ  ,  is  fixed  at  some 
multiple  Q  of  kR  .  For  purposes  of  the  present  study, 

Q  is  fixed  at  1000. 

The  rate  coefficients  and  &3  are  initially 

taken  as  unity,  but  then  varied  as  discussed  in  the  results 
section. 
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